Hanle effect driven by weak-localization 
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The influence of weak localization on Hanle effect in a two-dimensional system with spin-split 
spectrum is considered. We show that weak localization drastically changes the dependence of sta- 
tionary spin polarization S on external magnetic field B. In particular, the non-analytic dependence 
of S on B is predicted for III-V-based quantum wells grown in [110] direction and for [100]-grown 
quantum wells having equal strengths of Dresselhaus and Bychkov-Rashba spin-orbit coupling. It 
is shown that in weakly localized regime the components of S are discontinuous at B — 0. At low 
B, the magnetic field-induced rotation of the stationary polarization is determined by quantum 
interference effects. This implies that the Hanle effect in such systems is totally driven by weak 
localization. 



PACS numbers: 71.70Ej, 72.25.Dc, 73.23.-b, 73.63.-b 
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In recent years, the spin-related phenomena in semi- 
conductor nanostructures have been subject to intense 
study both in terms of fundamental physics and in view 
of applications in the field of spintronics Q]. The ul- 
timate goal of spintronics is to develop novel electronic 
devices that exploit the spin degree of freedom. The ef- 
fective manipulation of the spin in such devices requires 
that the characteristic spin lifetime be long compared to 
the device operation time. This is a challenging problem, 
especially for III-V-based semiconductor nanostructures 
where spin polarization relaxes rapidly due to Dyakonov- 
Perel spin relaxation mechanism jj. This mechanism is 
based on the classical picture of the angular spin diffu- 
sion in random magnetic field induced by spin-orbit cou- 
pling. In two-dimensional (2D) systems, the correspond- 
ing spin-relaxation time t$ is inversely proportional to 
the momentum relaxation time r 3]. As a consequence, 
in high-mobility structures which are most promising for 
device applications, ts is especially short. However, in 
some special cases, the relaxation of one of the spin com- 
ponents can be rather slow even in a system with high 
mobility. In particular, a number of recent researches 
El S El B 01 are devoted to GaAs symmetric quan- 
tum wells (QW) grown in [110] direction. In such wells, 
the random field is perpendicular to the QW plane and 
the normal to the plane component of the spin does not 
relax |3j or, more precisely, relaxes very slowly. Also, 
the random magnetic field might be parallel to a fixed 
axis in an asymmetric [100]-grown QW due to the in- 
terplay between structural and bulk spin-orbit coupling 
|lfl llll Il2l Il3l |l4j . Since in both cases one component 
of the spin relaxes slowly, these structures are especially 
attractive for spintronics applications. 

In this Letter, we discuss the dependence of the spin 
polarization in such structures on external magnetic field 
B. We assume that the spin is injected into the system 
with a constant rate, for example, by optical excitation 
[lij ]. The stationary spin polarization S is proportional 
to the product of injection intensity and the spin relax- 



ation time. The Hanle effect is that the external field 
modifies the stationary polarization. In particular, S(B) 
deviates from S(0) by an angle 9 which depends on the 
relation between spin precession frequency il = g/isB/fi, 
and the spin relaxation rate (here \xb = eh/2moc is the 
Bohr magneton, mo is the free electron mass and g is 
Lande g factor). We show that at low temperatures 9 is 
very sensitive to weak localization (WL) effects. Usually 
such effects are discussed in context of quantum correc- 
tions to the conductivity 0, 0] . Though the WL cor- 
rection is small compared to the classical conductivity it 
has attracted much attention due to its fundamental na- 
ture and anomalous behavior with external parameters 
such as magnetic field. Remarkably, the influence of WL 
on Hanle effect can not be considered as a small correc- 
tion to the classical result. We demonstrate that at low 
temperatures WL gives rise to a discontinuity in the de- 
pendence of 9 on B. As a result, at low B the Hanle effect 
is totally driven by WL. 

Physically, the non-analytic dependence of 9 on B is 
related to memory effects specific for WL. It is known 
that the WL is caused by interference of electron waves 
propagating along a closed loop in the opposite direc- 
tions [lij . Such interference process can be considered as 
a coherent scattering (additional to the Born scattering) 
involving a large number of impurities 03 (see Fig.©. 
The probability of the coherent scattering is proportional 
to the probability of the diffusive return l/Dt, where t 
is the time of the electron passage along the loop and D 
is the diffusion coefficient. In the absence of magnetic 
field, such scattering does not change the direction of 
the spin |l°J ]. Thus, electrons keep memory about ini- 
tial spin polarization during the time much larger than 
ts and the long- living tail 1/t in the spin polarization 
appears 0, |2(J • When the external magnetic field is ap- 
plied, the electron spin rotates with a frequency u) p + 17, 
where o; p is the momentum-dependent frequency of the 
spin precession in the spin-orbit-induced magnetic field. 
In the special case under discussion, u) p is parallel to a 



2 



fixed axis (say z-axis) for any p. Let us consider the sim- 
plest case, when fl also lies along z axis. In this case, spin 
rotation matrices, describing rotation of electron spin on 
the different segments of the closed loop commutes with 
each other and the spin rotation angles for clockwise and 
counterclockwise propagating waves are simply given by 
ip = Sit + J Q ujpdt' , tp' = Qt + J Q tOpidt' . The initial 
electron spinor \xo) is transformed to \x) = e~ la ' lp ^ 2 \Xo) 
and \x') = e~ l(Jzip / 2 |xo) f° r clockwise and counterclock- 
wise paths, respectively. Here a z is the Pauli matrix. 
In 2D systems, u> p is linear in p. Since for any closed 
path f*p(t')dt' = 0, we find tp = tp' = fit. The in- 
terference contribution to the electron spin-density ma- 
trix after a coherent scattering is given by |x)(x'l — 
e -ia z nt/2^ o ^^ e irr z Qt/2 _ jj e nce, the electron spin after 

the coherent scattering is rotated by an angle fit with re- 
spect to the spin before scattering. The deviation of the 
spin from the initial direction is proportional to sin(fii). 
The integration over t weighted with the probability of 
the coherent scattering yields 9 ~ J dtsm(flt)/t ~ fi/|fi|. 
In other words, the rotation frequency fi should be mul- 
tiplied by the effective rotation time l/|fi|, which is very 
long for small fi. Rigorous calculations (see below) give 
an additional factor A// in this result, which reflects the 
quantum nature of the phenomena: 9 ~ Afi/7|fi| (here A 
is the electron wave length and I is the mean free path). 
In contrast to coherent scattering, the classical rotation 
of spin is limited by the spin relaxation time rg, so that 
the classical contribution is 9 ~ Sirs- For firs <C A/Z, 
this contribution can be neglected and the Hanle effect 
is totally driven by WL. The discontinuity at fi = is 
smeared by inelastic scattering, which destroys phase co- 
herence between clockwise and counterclockwise propa- 
gating paths, thus limiting a time of the coherent spin 
rotation: t < r v , where is the phase breaking time. It 
worth noting that the above considerations do not work 
if w p changes its direction with a change of p. The reason 
is that the rotation matrixes corresponding to the differ- 
ent segments of the loop do not commute in this case. 
As a result, tp and ip' no longer equal to each other and 
depend not only on the total time spent on the loop but 
also on the positions of impurities 1, 2, ■ ■ • ,N. In such a 
situation, the external magnetic field leads to the spin 
decay rather than to the spin rotation [Igj . 

Next we develop a rigorous theory of Hanlc effect in 
the weakly localized regime. The Hamiltonian of a 2D 
system with a spin-split spectrum is given by 

H=^- + ^{u p + n]a + U(r). (1) 

Here p is the in-plane electron momentum, m is the elec- 
tron effective mass, a is a vector consisting of Pauli ma- 
trices, and U (r) is the impurity potential, which we as- 
sume to be short-ranged U(r)U(r') — j5(r — r') (here 
averaging is taken over impurity positions and the coef- 
ficient 7 is related to the transport scattering time by 




FIG. 1: A coherent scattering involving JV impurities. At 
point 1 the electron wave splits into two parts propagat- 
ing around the closed loop clockwise and counterclockwise. 
The initial electron spinor |xo) is transformed to \x) = 

e~ i& ^ /2 \xo) and | X ') = e _i *^' /2 |x) for clockwise and coun- 
terclockwise paths, respectively. 

t = h 3 /mj). The spin-orbit interaction is described by 
the term hu) p &/2. It can be separated into two parts 
(w p = up + up) related to so-called Dresselhaus jU 
and Bychkov-Rashba [2^] contributions. The Bychkov- 
Rashba coupling depends on the asymmetry of the QW 
confining potential. Its strength can be tuned by varying 
the gate voltage (2j|. The Dresselhaus term is present 
in semiconductors with no bulk inversion symmetry. In 
2D case the resulting spectrum splitting is linear in the 
electron momentum |3j, |22| . We assume that the random 
magnetic field is directed along z-axis 

u> p = (pa)z. (2) 

Here z is the unit vector along the direction of the random 
field and a is a constant in-plain vector. Equation 
implies that the spectrum splitting depends only on one 
component of momentum p a . This happens in symmetric 
[110]-grown QWs. In this case, z is normal to the QW 
plane |3|. For asymmetric [100] wells Eq. (0) can also take 
place if Bychkov-Rashba and Dresselhaus couplings have 
equal strengths O IU E EE 13 • For such wells vector 
z lies in the QW plane. The temperature is assumed to 
be low, so that t v be large compared to t$- 

To describe spin dynamics in the weakly localized 
regime we use the kinetic equation |l9l| . If the spin polar- 
ization is uniform in space, this equation looks as follows 

— = [wp + $7] x s + SJs + I, (3) 

at t 

where s = s(p) is the spin density in the momentum 
space, related to the averaged spin by S = J sd 2 p/ (2tt) 2 , 
s a = s — (s) is the anisotropic part of the spin density, 
Si = (s) is its isotropic part (here (...) stand for averag- 
ing over directions of the electron momentum), I is the 
constant source, which we assume to be monoenergetic, 
/ ~ S(E — Eq), and 8 J is the WL-induced correction to 
the Boltzmann collision integral [HI23. For E = E Q , 
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this correction reads (see Eq.(24) in Ref. (lif') 
XI 



SJs 



dt'W(t-t')s a (p,t'). (4) 



Here A = 2nh/^/2mE , I = Ty/2E /m and W(t) is the 
time- nonlocal scattering kernel j^. Using Eqs. © and 
(QJ we find a closed equation for the isotropic spin density 
in the stationary case 



f (1 - A)si - n x Si = I. 



(5) 



Here f = f s 1 is the spin relaxation tensor (tensor of 
inverse relaxation times) given by |3j 



E t fc = [8ik{w p ) - (wpiuipk)} t, 
and the matrix 

XI 



A 



dt W(t) 



(6) 



(7) 



describes WL correction (A^ <C 1) [26j. As follows from 
Eqs. J2J an< i ©? tensor T has two nonzero components: 
F xx = T yy = T = mE a 2 T .31]. We will assume that 
there exists a slow spin relaxation of the z component 
of the spin T zz — e <C T. Such relaxation arises due to 
the cubic in p terms, neglected in the Hamiltonian 0J, 
or due to other mechanisms of the spin relaxation. In 
the absence of the external magnetic field the scatter- 
ing kernel is given by Wik — Sik/AirDt 0. Therefore, 
Aife = (X/2ir 2 l)5ik ln(T v /V), where phase breaking time 
T(n is taken as an upper cut-off of the integral in Eq. J7J . 

Next we consider the dependence of A on the magnetic 
field for the case when B is parallel to z-axis. We start 
with discussing of the QW grown in [100] direction. Since 
in this case z lies in the QW plane, the external field does 
not affect the orbital motion of the electrons and the 
Zeeman term Mla z /2 commutes with the Hamiltonian. 
As a consequence, the solution of Eq. © with 1 = 
and B ^ can be obtained from solution at zero field: 
ss(i) = T(ilt)sB=o(t), where T is 3x 3 matrix, describing 
rotation around z-axis by the angle Qt. In order that 
ss{t) obeys Eq. ©, the scattering kernel has to be as 
follows: 



W(t) = 



T(Qt) 
A-nDt ' 



(8) 



Equations ©,0 and (JSJ) allows us to find the isotropic 
spin polarization. If I is perpendicular to the z-axis, we 
get 

S-(l + A c )S + (A s + O/r)zxS . (9) 
Here S = / Id 2 p/(27r) 2 r, S = / Sid 2 p/(27r) 2 , and 
A f 00 dt 



A r + lA s 
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(10) 



For Q/T <C 1, the angle between S and So is 

*(n)«S + Ai£r. (ID 



In this equation, f2/r stands for the classical contribu- 
tion, while the second term is due to WL effect. At small 
fields the classical contribution can be neglected. We see 
that WL contribution to 9 is a non-analytic function of 
the magnetic field. At zero field the function 9(0.) has a 
discontinuity 



(+0) - 0(-O) = X/2nl. 



(12) 



In deriving Eq. 1)11(1 . we neglected inelastic scattering. 
Such scattering destroys phase coherence, thus suppress- 
ing WL contribution. In order to take it into account, 
we introduce a factor exp(— t/r v ) in the integrand in 
Eq. (EH) - As a result, we find 



W) « f 



A 



27r 2 Z 



arctan(rir ¥ ,). 



(13) 



This equation shows that for low temperatures, when 
1/t v -c rA/2, the Hanle effect is still totally driven by 
WL at small fi. With increasing the temperature, the 
WL-induced discontinuity is smeared out. 

Let us now consider the system with symmetric QW 
grown in [110] direction, assuming again that O || z and 
I 1 z. The above considerations can be applied with a 
slight modification of Eq. (|10fl . Since in this case z is 
normal to the QW plane, the external magnetic field 
also effects orbital motion of the electrons leading to 
the suppression of WL 17]. Including into the integral 
in Eq. (|10fl an additional factor yt/ sinh(7i) j^], which 
accounts for destruction of WL by external field (here 
7 = 2eB D I he) , we obtain 



Q, AO f gX m 

f + 4r^M ta I "ST ^ 



(14) 



Next we find the stationary polarization for arbitrary 
directions of B and I. We choose x, y axes in such a way 
that fi x = 0. The classical solution (A = 0) reads 



i -o z /r n v /e 

n z /T l + Q 2 y /eT n z n y /eT 
-fty/e Vt z Vt y /eT T/e 



r Sp 



(15) 



(here we assumed ft -C T). To find A we write [13 

Wij(t) = ]T ( 1 \a l \(3)W Piae (0,t)(a\a J \9)/2, (16) 



where function W^ aS (r,t) obeys [H H |H IH 



K = ^ - D (V - 2ieA/hc + imZl&P + a (2) ]/2) 2 



(17) 
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Here <r$, ge , = (^|<x|/3')<W, trf^ gel = 5pp>(e\&\6') and 
A is the vector potential of external field. Using Eqs. JJJ, 
ijTfljl and fT7|l one can find 



A = 



A 



2tt 2 Z 



dt 



7 



sinh(7't) 



T(Q z t). 



(18) 



Here 7' = 2eHnD/hc and n is the unit vector normal 
to the QW plane. It is worth noting that the rotation 
matrix, entering Eq. 1)18(1 . depends only on z-component 
of the external field. Accounting for WL leads to the 
following replacement in Eq. 1)15(1 : 



r^r(i-A c ), n z 



ta. 



e^e(l-A ), (19) 



where A c = (A/2tt 2 Z) ln(t c /r), A = (A/2tt 2 Z) ln(i /r) 
and A s — (X/inl) ta,nh(\gmn z /8lmo\fln\) (here i^ 1 ~ 
max[\'Y\,T- 1 ], t- 1 ~ max[| 7 '|, \Vl z \]. Eqs. © and O 
are valid both for [110] and [100] orientations provided 
that fl z 3> AS2 Z ~ max[e, tly/T, 1/t v ]. One can show 
that the discontinuity of A s (at the point SI = ) is 
smeared over the frequency interval ACl z . 

A case when I is parallel to y axis is of particular in- 
terest. For Sl y ^ e and Q z <C TA S we find 



arctan 



02 



(20) 



We see that 6 > A s due to the factor rS\/(re+S!. 2 ) > 1. 
For rfl y /(Te + fly) ^> l/X the z-component of the spin 
will be larger than y-component. Hence, the quantum 
effects might lead to rotation of spin by a large angle. 

In the above calculations we assumed that I is homoge- 
nous. For slowly varying I, the derived equations relate 
S(r) with I(r) provided that the spatial scale of inhomo- 
geneity L is large compared to y/D/T = 1/ ma 30] . One 
can show that in the opposite case 1<L< 1/ma, these 
equations also valid relating Jc?rS(r) with J drl(r). 

Finally, we briefly discuss a possible experimental re- 
alization of the predicted effect. One of the most ef- 
ficient ways of the spin injection is the optical excita- 
tion of interband transitions with circular polarized light 
(lif . The observation of the weakly localized regime re- 
quires that t v 3> Ts- In optical experiments, the phase 
breaking is due to both inelastic scattering (caused by 
electron-phonon and electron-electron interactions) and 
recombination of electrons with holes. The recombina- 
tion of spin-polarized electrons with holes is suppressed 
in a n-type highly doped QW excited by low-intensity 
light. In this case, the number of spin polarized elec- 
trons is small compared to their total number and the 
holes are more probably to recombine with unpolarizcd 
electrons. Therefore, the stationary amount of holes is 
small and the recombination of polarized electrons with 
the holes can be neglected. The characteristic time of 
inelastic scattering will be large at low temperatures if 
the electrons are excited close to the Fermi level. In such 



a situation, the emission of the optical phonons is forbid- 
den and the phase breaking is due to electron-electron 
collisions. Since the rate of such collisions decreases with 
approaching to the Fermi level, r v can be tuned to be 
much longer than ts- 

To conclude, the theory of Hanle effect in a 2D system 
is developed for the weakly localized regime. At low ex- 
ternal magnetic fields the Hanle effect is totally driven by 
quantum interference effects. In the absence of inelastic 
scattering, the components of the spin polarizations are 
discontinuous as functions of the external field. 

We are grateful to K. V. Kavokin for useful discussions. 
This work has been supported by RFBR, a grant of the 
RAS, a grant of the Russian Scientific School 2192.2003.2, 
and a grant of the foundation "Dinasty"-ICFPM. 



[1 

[2 

[3 

[4] 
[5] 
[6] 

ir 

[8 
[9 
[10 

[11 

[12 

[13 
[14; 
[IS 
[16 

[i7; 
[is; 

[19 

[20 

[21 

[22 
[23 
[24 

[25 
[26 



Semiconductor Spintronics and Quantum Computa- 
tion, eds. D.D. Awschalom, D. Loss, and N. Samarth 
(Springer- Verlag, Berlin, 2002). 

M.I. Dyakonov and V.I. Perel', Sov. Phys. Solid State 13, 
3023 (1972). 

M.I. Dyakonov and V.Yu. Kachorovskii, Sov. Phys. Semi- 
cond. 20, 110 (1986). 

Y. Ohno et al, Phys. Rev. Lett. 83, 4196 (1999). 

T. Adachi et al, Physica E (Amsterdam) 10, 36 (2001). 

W.H. Lau and M.E. Flatte, J. Appl. Phys. 91, 8682 

(2002) . 

O.Z. Karimov et al, Phys. Rev. Lett. 91, 246601 (2003). 

K.C. Hall et al., Phys. Rev. B 68, 115311 (2003). 

S. Dohrmann et al, Phys. Rev. Lett. 93, 147405 (2004). 

F. G. Pikus and G.E. Pikus, Phys. Rev. B 51, 16928 
(1995) 

N.S. Averkiev and L.E. Golub, Phys. Rev. B 60, 15582 
(1999) 

A. A. Kisilev and K.W. Kim, Phys. Status Solidi (b) 221, 
491 (2000) 

J. Schliemann, J.C. Egues, and D. Loss, Phys. Rev. Lett. 
90, 146801 (2003). 

J. Schliemann and D. Loss, Phys. Rev. B 68, 165311 

(2003) . 

Optical orientation, eds. F. Meier and BP. Za- 
kharchenya, (North Holland, Amsterdam, 1984). 
LP. Gor'kov, A.I. Larkin, and D.E. Khmelnitskii, JETP 
Lett. 30, 248 (1979). 

For review, see P.A. Lee and T.V. Ramakrishnan, Rev. 
Mod. Phys. 57, 287 (1985). 

A. P. Dmitriev, I.V. Gornyi, and V.Yu. Kachorovskii, 
Phys. Rev. B 56, 9910 (1997) 

I.S. Lyubinskiy, V.Yu. Kachorovskii, Phys. Rev. B, in 
press 70, N. 16 (2004), cond-mat/0404391 
A.G. Mal'shukov, K.A. Chao and M. Willander, Phys. 
Rev. B 52, 5233 (1995). 

G. Dresselhaus, Phys. Rev. 100, 580 (1955). 

Yu. Bychkov and E. Rashba, JETP Lett. 39, 78 (1984). 
J. Nitta et al., Phys. Rev. Lett. 78, 1335 (1997). 
Assuming that the cyclotron frequency is small compared 
to 1 /t, we neglect the term, which describes the momen- 
tum rotation by external field. 

To obtain Eq. Q) from Eq. (24) of Ref. [r| we used the 

following property s a (— p,i) = — s a (p,t). 

In deriving of Eq. (J^J, we took into account that W is 



5 



proportional to the rotation matrix (see Eqs. JHJ, USt - ). 
[27] S. Chakravarty and A. Schmid, Physics Reports, 140, 
195 (1986) 

[28] W. Knap et al, Phys. Rev. B 53 3912 (1996) 

[29] S.V. Iordanskii, Yu.B. Lyanda-Geller, and G.E. Pikus, 



JETP Letters 60, 206 (1994). 
[30] Note that 1/ma is a characteristic scale of non-decaying 
spin oscillations which might exist in the system under 
discussion ( see Ref. Hal ). 



